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We use Time Resolved Correlation (TRC), a recently introduced light scattering method, 
to study the dynamics of a variety of jammed, or glassy, soft materials. The output of 
a TRC experiment is ci(t,r), the time series of the degree of correlation between the 
speckle patterns generated by the light scattered at time t and t + r. We characterize the 
fluctuations of cj by calculating their Probability Density Function, their variance as a 
function of the lag t, and their time autocorrelation function. The comparison between 
these quantities for a Brownian sample and for jammed materials indicate unambiguously 
that the slow dynamics measured in soft glasses is temporally heterogeneous. The analo- 
gies with recent experimental, numerical and theoretical work on temporal heterogeneity 
in the glassy dynamics are briefly discussed. 

Keywords: Light scattering, fluctuations, slow dynamics, glass, jamming, dynamical 
heterogeneity. 

1. Introduction 

Dynamic light scattering (DLS), also known as photon correlation spectroscopy, is 
a popular technique that provides information on the dynamics of soft condensed 
matter systems [1]. It is widely used in physics, chemistry, and biology to investigate 
systems as diverse as colloidal particle suspensions, emulsions, polymer or protein 
solutions, surfactant phases, and gels. In the past years, new methods and tech- 
niques have been developed that have greatly extended the range of applicability of 
dynamic light scattering. The advent of coherent, high brilliance X-ray sources has 
allowed photon correlation spectroscopy experiments to be performed using X-ray 
radiation [2], instead of the laser beam used for DLS. Since in a scattering experi- 
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ment the wavelength of the radiation sets the length scale over which the dynamics 
is probed, X-ray photon correlation spectroscopy (XPCS) is a significant extension 
of DLS. The interpretation of both DLS and XPCS data is based on the assump- 
tion that the scattering is sufficiently weak, so that virtually all photons collected 
by the detector are singly scattered and the distortions induced by multiple scat- 
tering can be neglected. Diffusing Wave Spectroscopy (DWS) [3], on the contrary, 
allows one to study very turbid samples, for which multiple scattering dominates, 
thus extending the range of samples to which dynamic scattering methods can be 
applied. 

Quite generally, in a scattering experiment the information on the dynamics is 
extracted from the time fluctuations of the scattered intensity, I(t). The inten- 
sity is usually measured by a point detector, such as a phototube; its fluctuations 
are quantified by calculating the time autocorrelation function of I(t), defined by 
C?2(t) = (I(t)I(t+ T ))t, where (...)t denotes an average over time 1 . The time average 
is necessary to obtain sufficiently good statistics: typically, data have to be collected 
over a period as long as 10000 times the relaxation time of G2(t). Because of the 
extended time averaging, traditional scattering methods can only provide informa- 
tion on the average dynamics of a system, with no time resolution. The lack of time 
resolution, together with the need of collecting data over a long period, makes tradi- 
tional scattering techniques difficult to apply to glassy — or jammed [4] — systems, 
for which there is currently a great interest (for a review, see for example Ref. 5). 
In fact, for these systems the relaxation time can be as long as several thousand 
seconds, thus making time averaging practically impossible. Moreover, their dy- 
namics is often non-stationary; therefore, time-averaged correlation functions do 
not allow the evolution of the dynamics to be followed precisely. The introduction 
of the so-called multispeckle technique has addressed these issues. In a multispeckle 
experiment, the point detector is replaced by the 2-D array of pixels of a CCD cam- 
era [6-8]. Each CCD pixel acts as an independent detector, for which an intensity 
time autocorrelation function is calculated. In multispeckle experiments, the time 
average that is usually performed in traditional light scattering is replaced by an 
average over the CCD detector: G 2 (i, i") = (I p (t)I p (t + t)) p , where the subscript p 
refers to the CCD pixels. In most works, G%(t, r) is also averaged over a short time 
window [8] or over several repetitions of the experiment [9] , in order to reduce the 
measurement "noise" . 

Recently, we have proposed a novel scattering method, termed Time Resolved 
Correlation (TRC), where valuable physical information on the dynamics are ex- 
tracted precisely from this "noise", i.e. from the time fluctuations of G^^r) [10]. 
In particular, we have shown that, contrary to traditional scattering methods, TRC 
allows temporally heterogeneous dynamics to be discriminated from homogeneous 
dynamics 2 . Because there is growing evidence that dynamical heterogeneity is a 
key feature of the slow dynamics of glassy systems, TRC appears as a very promis- 

1 Although in principle information on the dynamics may also be extracted from the spectral 
properties of the scattered light, the analysis in the temporal domain has become the standard 
procedure. 

2 As it is common practice in the literature on glassy systems, in this paper we use the adjective 
"heterogeneous" to indicate the presence of fluctuations in the dynamics other than the inevitable 
noise associated with the measuring process. Note that, with this definition, both stationary and 
non-stationary dynamical processes may be temporally heterogeneous. 
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ing tool for investigating the slow relaxation of these systems. In this paper, we 
review the principles of TRC and we propose some ways of analyzing TRC data. 
The paper is organized as follows: Sect. briefly reviews the basic features of a 
scattering experiments and the way TRC data are collected. In Sect.|3|we develop 
the statistical analysis used to characterize the fluctuations of the TRC signal and 
we apply it to data collected for systems whose dynamics is homogeneous. In the 
following Section, the same tools are used for systems that exhibit temporal het- 
erogeneities. Section concludes the paper by summarizing the main results and 
briefly compares them to other experimental and numerical works on systems whose 
dynamics is temporally heterogeneous. 

2. Time Resolved Correlation 

The TRC method can be applied both to the single scattering geometry and to DWS 
experiments. In all cases the multispcckle detection is adopted: a CCD camera is 
used to take a time series of images of the speckle pattern, the characteristic grainy 
picture generated by the interference of the photons scattered by the sample. 

Typical scattering setups are sketched in Fig. for single (a) and multiple 
(b) scattering measurements; XPCS instruments are similar, the laser source being 
replaced by synchrotron radiation. A detailed discussion on how the speckle pattern 
fluctuates because of the sample dynamics is beyond the scope of this paper and 
can be found in Ref. 1,3. For our purposes, we just note that any change in the 
position of the scatterers modifies the relative phase of the scattered photons and 
thus the intensity pattern of the speckles generated by their interference. Using the 
notation of Ref. 10, we quantify the change in the speckle pattern by introducing 
c/(t, r), the normalized intensity correlation between speckle images at time t and 
t + t: 

r(tr )- G ^ T ) 1- (Ip(t)I p (t + T)) p 

n ' J (Ip(t)) P (I P (t + T)) p (I p (t)) p (I p (t + r)) p ' {l 

where the average is over CCD pixels. The normalization in Eq. guarantees that 
cj (i, r) is not affected by any change in the incoming beam intensity. Because of 
the statistical properties of the speckle pattern, < cj(t,r) < (3, where f3 < 1 is a 
prefactor that depends on the setup geometry [1, 11]. 

Figure [5] shows a 2-D grey scale representation of cj as a function of both time 
t and time lag r (data were taken in the DWS geometry for a colloidal paste [12]). 
White corresponds to the highest degree of correlation, cj = (3, while black corre- 
sponds to uncorrelated pairs of speckle patterns, cj — 0. As general trend, for any 
given t ci decreases with increasing time lag r, since the longer the lag, the greater 
the change of the system configuration and thus that of the speckle images. This is 
exemplified by the right panel that shows cj(t = 3220 sec, r), obtained by cutting 
the 2-D data along the vertical dashed line. Note that the time average of c/(t, r) 
yields the normalized intensity autocorrelation function, 52 (t) — 1 , usually measured 
in light scattering experiments 3 . The time fluctuations of cj measured at a constant 
lag t = 1360 sec, obtained from a cut of the 2-D data along the horizontal solid line, 
are shown in the top panel. Sharp drops of correlation that depart from the typical 



3 In traditional light scattering measurements, the average is exclusively over time when using a 
point detector, over both pixels and time when using a CCD camera. 
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Fig 1. Schematic representation of a) a single scattering apparatus, b) a multiple scattering (DWS) 
setup in the transmission geometry. In a), 9 is the scattering angle, selected by the diaphragm 
in front of the CCD . The lens in b) conjugates the plane indicated by the dashed line with the 
CCD chip. In both cases, the CCD camera is used to record the speckle pattern generated by the 
interference of the scattered photons. 
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Fig 2. Bottom left panel: 2-D grey scale representation of Cf(t,r), the degree of correlation 
between speckle patterns recorded at time t and t + r. White is the highest correlation, black 
the lowest. The sample is a concentrated colloidal paste. Right panel: two-time intensity auto- 
correlation function obtained by cutting cj along the vertical dashed line. Top panel: temporal 
fluctuations of cj at a fixed lag time, obtained from a cut along the horizontal solid line. Note the 
sharp drops of correlation, indicative of sudden, temporally localized rearrangements. 
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fluctuations of the signal are clearly visible. These drops suggest that occasionally 
the change in the speckle pattern, and thus in sample configuration, is much larger 
than average: they are therefore indicative of a temporally heterogeneous dynamics. 
In the rest of the paper, we will focus on the statistical analysis of the fluctuations 
of Ci(t, t) at fixed lag r, first for samples that exhibit homogeneous dynamics, and 
then for systems whose dynamics is temporally heterogeneous. 

3. Statistical Analysis of TRC Data for temporally homogeneous dy- 
namics 

In order to better appreciate the signature of heterogeneous dynamics in TRC 
data, we first characterize the statistics of the fluctuations of cj(t,r) in the case 
of temporally homogeneous dynamics. The data presented in this Section are taken 
in the single scattering geometry, at a scattering angle 9 — 90 deg 4 . The sample is a 
diluted suspensions of particles undergoing Brownian motion: polystyrene spheres 
of radius r = 530 nm were suspended at a volume fraction <f> — 3.7 x 10 -5 in 
almost pure glycerol (approximately 99.9/0.001 w/w glycerol/ water). The sample 
temperature was controlled and fixed at T = 15 ± 0.1 °C. 

Figure[3]a) shows the time-averaged intensity autocorrelation function 52 (t) — 1 
obtained from the TRC data for the Brownian particles. The data are well fitted 
by an exponential decay /3exp(— t/t s ) (line in the figure), as expected for a dilute 
suspension of monodisperse spheres [1], with t s = 2.35 sec. The right panel shows 
the corresponding ci(t,r) as a function of time t for four time lags, ranging from 
t/t,s = 0.21 to t/t s = 6.38. For Brownian particles, we expect the dynamics to 
be temporally homogeneous. Accordingly, the loss of correlation between pairs of 
speckle images separated by a given time lag r should not depend on the time the 
first image of the pair is taken. In other words, for a fixed lag r, ci(t,r) should be 
constant. This is indeed the case, except for the small fluctuations visible in in the 
right panel, which we attribute to the CCD noise and the finite number of pixels, 
and therefore independent speckles, over which c/ is averaged. In the following we 
shall refer to both sources of noise as to the "measurement noise" . 

In order to verify this hypothesis, we analyze in some detail the temporal fluc- 
tuations of c/. Figure 0] shows the Probability Density Function (PDF) of c/ for 
the same time lags as in Fig. [21 right panel. At all lags, the data are well fitted by 
a Gaussian distribution, as shown by the lines. This behavior is consistent with the 
hypothesis that the fluctuations are due to the measurement noise. Note that the 
width of the distributions decreases for increasing lag r, together with the mean 
value 52(7") — 1 = (cj)t- To better investigate the lag dependence of the amplitude 
of the fluctuations, we plot in the left panel of Fig.[5]the variance of c/, var(c/), as 
a function of r normalized by the mean relaxation time t s . At the smaller lags, the 
decay of the amplitude of the fluctuations closely follows that of 52 (t) — 1, whereas 
for t > t s var(cj) saturates to a plateau. This behavior suggests the following form 
for the variance of Cj; var(cj) = A + B{{a) t + l) 2 , where the second term of the 
r.h.s. represents a noise contribution proportional to Gi(t, t) / {I p {t)) p {I p (t + r)) p , 
the signal that is measured in a TRC experiment. The first term is the noise of 
C/ when all correlation between speckle images is lost, i.e. the fluctuations of the 



4 Similar results are obtained in DWS 
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Fig 3. For a dilute suspension of Brownian particles arc shown: left panel: the intensity auto- 
correlation function <?2(t) — 1; right panel: cj(t,r) as a function of t for four time lags r. The 
correlation function shown in the left panel was obtained by time-averaging the TRC data plotted 
in the right panel (for clarity, not all the time lags are shown in the right panel). 




Fig 4. Probability Density Function of the temporal fluctuations of cj(t, r), for the data shown in 
Fig. EJb) (symbols: from left to right, r = 0.5, 1.5, 3, and 15 sec). The lines are Gaussian fits to 
the data. Note that width oh the PDF's — and thus the standard deviation of the fluctuations — 
decreases at larger lags, together with the mean of c/. 



A. Duri, P. Ballesta, L. Cipelletti, H. Bissig and V. Trappe 



base line for t — » oo. As it will be discussed in more detail in reference to Fig. 
this contribution is mainly due to the finite number of pixels. The right panel of 
Fig. shows var(cr) vs ((c/) t + l) 2 (symbols). The data are in good agreement 
with the linear behavior var(cj) = A + B((cj) t + l) 2 introduced above, except for 
the downturn visible at the smallest values of (c/)t (the straight line is a linear fit 
for ((c/)t + l) 2 > 1.3). The origin of this deviation most likely lies in correlations 
between G%(t, r) and (I p (t)) p . We are currently investigating the origin of this issue. 

The measurement noise is expected to increase as the number of pixels over 
which the intensity correlation is averaged, n p i X , decreases. More precisely, the 
variance of the CCD electronic (or dark) noise and read-out noise should scale as 
n pix > because of the Central Limit Theorem. Similarly, even in the absence of CCD 
noise, all pixel-averaged quantities in Eq. JJ) will fluctuate as the speckle pattern 
is renewed over time, because n P i X is finite [11]. By invoking again the Central 
Limit Theorem, we expect the variance of these fluctuations to exhibit the same 
n P ix scaling as those due to the CCD noise. Therefore, the overall variance of the 
measurement noise should decay with n P i X as a power law with exponent -1. The 
left panel of Fig.[|)]shows that indeed this behavior holds for all r, both smaller and 
larger than the mean relaxation time r s , thus further supporting the hypothesis 
that the fluctuations of cj for temporally homogeneous systems are due to the 
measurement noise. 

The PDF and the variance of cj that we have so far introduced are not sensitive 
to the detailed way c/ fluctuates over time, but only to its probability to assume a 
given value. To investigate the time behavior of the fluctuations, we calculate the 
time autocorrelation function of c/(i,r), defined as 

p„ ( s fA .s {(ci(t,T) ~ ( Cl ) t )( Cl (t+At,T) - ( Cl ) t ))t , . 

Corr (c/, T ) {At) = , 2 

{{Cl{t,T) - {Cl)tY)t 

where we have added the index r to designate the time lag for which cj is calculated. 
Note that with the above definition Corr(c/ iT )(0) = 1 and Corr(cj jT )(Ai) = for 
uncorrelated signals (e.g. for At — > oo). The quantity introduced in Eq. J5J is a 
fourth-order time correlation function, since cj itself is a time correlation function. 
In order to avoid any confusion, in the following we shall refer to the argument At 
of Corr(c/ iT ) as to the time delay, as opposed to the time lag r for cj. In the 
right panel of Fig. Corr(c/ !T )(Ai) is plotted as a function of At normalized by 
the relaxation time r s , for the same four time lags as in the previous figures of this 
section. As it can be seen, Corr(c/, r )(Ai) decays to on a time scale comparable to 
the relaxation time of the sample, r s . Indeed, the data can be fit by an exponential 
relaxation, whose characteristic time varies between 0.21r s and 0.31r s (for clarity, 
only the fit for the data collected at t = 1.5 sec is shown, as a solid line, in the 
figure). At first sight this behavior may seem in contrast with the assumption that 
here the fluctuations of c/ are due to the measurement noise, which is expected to 
be uncorrelated. Indeed, this is the case for the CCD noise, as we tested by taking 
a series of dark images. However, we recall that the fluctuations of the quantities 
(I p {t)I p {t + t)) p and (I p ) p used to compute cj are due not only to the CCD noise, 
but also to the finite number of speckles that are sampled when averaging over the 
CCD pixels. This latter contribution to the noise of cj(i, t) is correlated over a 
time scale comparable to r s , for all time lags r. To demonstrate it, let us consider 




Fig 5. Left: lag dependence of the variance of the temporal fluctuations of cj for a suspensions 
of Brownian particles. Right: same data plotted parametrically vs the squared mean of the TRC 
signal, ((cj)t + l) 2 . Except for the lowest values of ((cj)t + l) 2 , the data follow a linear behavior, 
as shown by the straight line. 
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Fig 6. For a suspension of Brownian particles are shown: Left panel: double logarithmic plot of 
the variance of the temporal fluctuations of c; , as a function of the number of processed pixels, 
and for various time lags. For all r, the data are well fitted by a power law decay with exponent 
—0.98 ± 0.01, very close to what is expected for the measurement noise because of the Central 
Limit Theorem. Right panel: time autocorrelation function of the fluctuations of cj, for various 
lags (symbols). The line is an exponential fit to the data for r = 1.5 sec (the fit for the other values 
of t are omitted for clarity). The correlation in the fluctuations stems from the finite number of 
pixels over which cj is averaged, as discussed in the text. 
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as an example the fluctuations of Ci(t,r) for r S> t s , so that (ci(t, r)) f = 0. A 
measurement of c/ at a given time t will yield a value e that is in general slightly 
different from zero, because c/ is averaged over a finite number of pixels. If the 
measurement is repeated at a later time t + At, with At ^ t s , the outcome will be 
essentially the same, since the system configuration is almost frozen on a time scale 
At <C t s and thus the speckle patterns recorded at time t + At and t + t + At are 
almost unchanged with respect to those recorded at time t and t + T, respectively. 
Therefore, cj(t + At,r) is highly correlated to cj(t, r). As the delay is increased, 
the evolution of the speckle images during At will become increasingly relevant and 
cj(t + At, t) and C[(t,r) will become less and less correlated. Since the speckle 
fields are renewed completely after a few t s , one expects Corr(c/ iT )(At) to decay 
with a characteristic time close to t s , as observed in the right panel of Fig. H3 (note 
that similar arguments can be invoked for all r). 

4. Statistical Analysis of TRC data for temporally heterogeneous dy- 
namics 

We now turn to TRC data collected for jammed systems and show that the statis- 
tical properties of c/ are very different from those reported in the previous section 
for dilute Brownian particles. Although wc will focus in turn on different systems, 
such as a micellar polycrystal, the colloidal paste cited above, and a shaving cream 
foam, we stress that most of the features that will be shown appear to be quite 
general. As we shall discuss in the following, these differences indicate that the 
slow dynamics of jammed systems is temporally heterogeneous. 

Figure[3shows TRC data collected in the single scattering geometry {9 = 45 deg) 
for the slow dynamics of the network of defects of a micellar polycrystal (more details 
on the sample and its average dynamics can be found in Ref. 13 and the references 
therein). For this system, the time-averaged intensity autocorrelation function is 
a "compressed" exponential: 32(1") — 1 = /3cxp[— (t/t s ) p ], with r s = 3160 sec and 
p w 1.7. Similarly to the case of the colloidal paste data shown in the upper panel of 
Fig. the fluctuations of cj are very different from those measured for a Brownian 
suspension (compare to Fig.|3 right panel). As discussed previously, the large drops 
of ci(t,r) below its mean value indicate that occasionally the sample undergoes a 
rearrangement that is much larger than the typical change in configuration over 
a time r. Therefore, the TRC data unambiguously show that the dynamics of 
both the colloidal paste and the micellar polycrystal is intermittent. The PDF 
of ci{t,r = 750 sec) is shown in the right panel. Because of the drops of the 
signal visible in the left panel, the PDF severely deviates from the Gaussian shape 
observed for the Brownian sample (compare to Fig. and is strongly skewed to 
the left, the left tail being close to an exponential decay. Interestingly, the shape of 
the PDF depends on the time lag at which cj is measured. Without discussing the 
details of this dependence, we remark that for r ^ r s the PDF recovers a Gaussian 
shape (data not shown). This behavior can be understood by assuming that the 
average time between rearrangements is much shorter than the mean relaxation 
time r s . Accordingly, a large number of rearrangements occur between any two 
images taken at a lag r ;§> t s and the PDF is Gaussian because of the Central Limit 
Theorem. In other words, on these time scales the dynamics is due to such a large 
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Fig 7. For a micellar polycrystal arc shown: Left panel: time variation of cj(t,r) for several lags 
(the curves are labelled by the normalized lag t/t s ). The relaxation time of g2 — 1 is t s = 3158.5 
sec. Right panel: for the same data, PDF of cj for t/t s = 0.237. Note the highly skewed shape 
of the PDF, due to the marked drops of correlation visible in the left panel. For the longest r, a 
Gaussian shape of the PDF is recovered (not shown in the figure). 
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number of individual rearrangement events that it can be regarded as continuous, 
much as for a Brownian suspension the dynamics is continuous because it arises 
from the motion of a very large number of particles, even though each individual 
particle moves in a jerky fashion. 

We continue the comparison between the statistical properties of TRC data 
for jammed systems and those for fluid systems by contrasting the t dependence 
of the variance of cj for the colloidal paste introduced in Sect. [21 to that for the 
diluted Brownian suspension, shown previously in Fig. [SJ As can be seen in Fig. EI 
for the colloidal paste var(c/) is peaked around r « t s (for this sample, 52 (t) — 
1 « /3exp(— t/t s ), with r s = 1580 sec), instead of decaying with r, as for the 
Brownian particles. We stress that the peaked shape is conserved even when plotting 
parametrically var(c/) vs ((ci)t + l) 2 (data not shown), contrary to the case of the 
Brownian suspension (see Fig. [SJ right panel). Moreover, the peak value of the 
variance is 2 orders of magnitude larger than the measurement noise at large r, 
whereas for the Brownian sample the difference between the smallest and largest 
value of var(c/) is only about a factor of 3. Therefore, both the presence of a 
peak and the larger amplitude of the fluctuations clearly differentiate the dynamics 
of the pastes from that of the Brownian particles, even though the shape of the 
time-averaged correlation function is the same for both systems, i.e. an exponential 
decay. The presence of a peak in var(cj) can be understood, at least qualitatively, 
by assuming that the characteristic time of the relaxation of cj fluctuates with 
time, r s = T s (t), rather than being constant. To illustrate schematically this point, 
we show in the right panel of Fig. [H] two exponential decays with slightly different 
relaxation times (in arbitrary units), t s = 1 and 1.5, respectively. These two curves 
are meant to be representative of the typical behavior of cj vs r measured at 
two different times t\ and £2- As indicated by the arrows, the resulting change in 
correlation at a given lag is maximum for r w t s and decreases both at small and 
large lags, in qualitative agreement with the curve shown in the left panel. For our 
jammed systems, the physical origin of the fluctuations of r s most likely lies in the 
fluctuation of the number, or the size, of the regions that rearranged by each event. 
More work is currently being carried on to explore quantitatively this concept. 

To investigate how the fourth-order correlation function Corr(cj. r )(Ai) behaves 
when the dynamics is temporally heterogeneous, we analyze TRC data collected 
in the DWS transmission geometry for a shaving cream foam. Previous work by 
Durian et al. has shown that the dynamics of a foam is intermittent, due to bubble 
rearrangements that are localized both in time and in space [14]. Therefore, the 
foam represents an ideal system for testing the TRC method. In our experiment, we 
find that the time-averaged correlation function #2 — 1 is well fitted by an exponential 
decay with characteristic time r s = 0.29 sec, in agreement with Ref. 14. The 
left panel of Fig. shows Corr(c/. r ) vs. the normalized delay At/r s for various 
lags r. In sharp contrast with the behavior of the fourth-order correlation for 
the Brownian sample (see Fig. |HJ right panel), for the foam Corr(c/. T ) strongly 
depends on the lag at which cj is calculated, the decay time of Corr(cr jT ) increasing 
with r. To characterize this dependence, for each r we extract from the initial 
relaxation of Corr(c/ iT ) a (normalized) characteristic time, T^, which we define 
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Fig 8. Left panel: r dependence of the variance of cj for a dense colloidal paste. Data are taken 
in the DWS transmission geometry. Note the peak of var (cj) for r « r s and the large value of the 
fluctuations at the peak, compared to that for t ^s> t s . Right panel: schematic representation of 
the impact of a fluctuation of the relaxation time r s on the value of cj. As shown by the arrows, 
the change of cj at a fixed lag is maximum for r f» t s = 1. 
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Fig 9. For DWS measurements on a shaving cream foam are shown: Left panel: the time 
autocorrelation of cj, for various lags t, as a function of the normalized delay At/r s . The curves 
are labelled by t/t s . Contrary to the case of Brownian particles, there is a strong r dependence of 
the decay rate of Corr(cj T )(At). Right panel: r dependence of the relaxation time of the curves 
shown in the left panel (symbols, see text for details). For the shortest lags, Fgo is proportional 
to r, as indicated by the straight line. 
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as the (normalized) delay for which Corr(c/ iT ) = C t h = 0.90 5 . We plot Tgo as a 
function of r in the right panel of Fig. El as it can be seen, for small lags Tgo is 
proportional to r, whereas it deviates from this linear trend as r approaches the 
relaxation time of 52 — 1- Therefore, for the foam the fluctuations of ci(t,r) are 
correlated over a time scale proportional to r, at least as long as the time lag is not 
too large. 

In order to explain this behavior, we represent schematically the evolution of 
a system whose dynamics is intermittent in Fig. The rectangles at different 
grey levels represent different configurations of the system: the evolution occurs 
via discrete "jumps" — or rearrangements — from a frozen configuration to the next 
one, rather than through a continuous change. The life time of the configurations 
is assumed to vary randomly. The degree of correlation between pairs of speckle 
images drops whenever the two images are taken for different configurations, while 
it stays high within the same configuration (see the sketch of c/ below the time 
axis in Fig.^J for two values of r). This implies that the dips of cj have a width 
equal to the time lag r, and that their raising edge is aligned with the change in 
configuration, thus allowing one to identify precisely the time at which the rear- 
rangement event occurred. Accordingly, for any given r, the TRC signal consist of 
a random sequence of downward pulses of width r. The autocorrelation of such a 
signal decays over a time scale of the order of the width of the pulses, i.e. over a 
time of order r. Therefore, this rather crude picture provides some explanation to 
the proportionality between the decay time of Corr(c/ jT ) and r that is observed, 
at small r, for the foam. Note that as the lag increases the likelihood of two or 
more rearrangements to occur during a time r becomes increasingly relevant. These 
multiple rearrangements would alter the shape of the dips of c/, thus reducing the 
decay time of Corr(cj iT ) with respect to what is expected by extrapolating the small 
r regime, as observed in Fig. |5] (right panel). Of course, a more realistic description 
of the dynamics should take into account the possibility that the rearrangement be 
more gradual and the dynamics within a given configuration not completely frozen. 
Moreover, in a macroscopic scattering volume many rearrangements are likely to 
occur at any given time, so that in practice r would be always longer than the mean 
life of the configuration. This being the case, a more appropriate approach would be 
to relate the fluctuations of ci(t, r) to those of the number of rearrangements during 
a time r, rather to the occurrence of single events. We are currently exploring this 
scenario. 

5. Conclusions 

We have shown that, contrary to traditional light scattering techniques, TRC ex- 
periments have the capability of discriminating between continuous and temporally 
heterogeneous dynamics. Among the numerous ways one could analyze the TRC 
time series, three methods seem particularly interesting: the characterization of 
the PDF of the fluctuations of cj, the study of the r dependence of the width of 
the distributions, var(cj. T ), and the investigation of the temporal correlations via 
Corr(c/ iT ). Other possible approaches may include a wavelet analysis of c/ and the 
investigation of the fluctuations in the framework of chaos theory. Quite generally, 



5 Similar results are obtained for different Cth> provided that C t h > 0.7. 
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Fig 10. Schematic representation of the evolution of a system, such a foam, whose dynamics 
is intermittent: each rectangle represents a (slightly) different configuration of the system. The 
time evolution of C/ for two values of the lag, n and T2, is shown below the time axis: cj drops 
whenever the two speckle images that are correlated are taken for different sample configurations. 
Note that the width of the dips is equal to the time lag. 
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TRC measurements indicate that the slow dynamics of soft jammed materials is 
temporally heterogeneous, in striking analogy with recent numerical and experi- 
mental work on hard condensed matter glasses (for a review, see Refs. 15, 16). We 
conclude by a brief, and partial, overview of some of the recent work related to the 
investigation presented in this paper, either because of the techniques that are used 
or in view of the results that are obtained. Dixon and Durian have introduced a 
method termed Speckle Visibility Spectroscopy (SVS), where the fast dynamics of a 
multiply scattering sample is studied by measuring with a CCD the time-dependent 
contrast of the speckle pattern, defined by (I P (t) 2 ) p / (I p (t))p — 1 [17]. Note that, in 
our notation, this quantity corresponds to ci(t, r = 0). If the speckles fluctuate on a 
time scale shorter than or comparable to the CCD exposure time, the image will be 
blurred and the contrast reduced. Therefore, SVS provides information on the dy- 
namics on a time scale comparable to the exposure time (typically from a fraction of 
msec to a fraction of sec), while in the TRC method the shortest time scale is set by 
the CCD acquisition rate (typically a few tens of images per second) and the longest 
one can be as large as several hours. Other similarities exist between the TRC and 
the work by Lemieux and Durian, who use fourth-order intensity correlation func- 
tions to measure intermittent dynamics [18,19]. Although the quantity Corr(c/ !r ) 
introduced here is also a fourth-order intensity correlation function, we note that in 
Refs. 18, 19 a point detector is used and no pixel average is performed, as opposed 
to the way cj is calculated. An extended time average is thus required to reduce 
the measurement noise; contrary to TRC, no time resolution can be obtained: it is 
possible to access some of the average statistical properties of the fluctuations, not 
their detailed temporal behavior. An interesting parallel can be made between the 
work presented here and the investigation of temporally heterogeneous dynamics in 
2-D electron systems [20,21]. Following a method originally designed by Weissman 
in the context of spin glasses [22] , a time series of the fluctuations of the resistance 
(or the conductivity) is divided into shorter segments. For each segment, the power 
spectrum of the fluctuations is calculated. Finally, the so-called second spectrum 
is obtained by calculating the power spectrum of the time series of spectra calcu- 
lated for the individual data segments. Because of the relationship between power 
spectrum and autocorrelation function of a signal, the fourth-order autocorrelation 
function Corr(c/ j7 -) calculated here for TRC represents the time-domain analogue 
of the analysis proposed in Ref. 22. The PDF of the fluctuations of the dynamics of 
a polymeric glass, as obtained from dielectric measurements, has been reported by 
Buisson et al. [23,24]. Skewed, non-Gaussian PDFs are observed, similar to that re- 
ported in Fig.0 Non-Gaussian distributions are also observed in simulations of spin 
glasses, where the dynamics was averaged over subregions instead of over the whole 
simulation box [25,26]. Interestingly, for both our data (Fig.EJ) and the simulations 
on spin glasses, the shape of the PDF is very close to the universal non-Gaussian 
PDF proposed for the fluctuations in highly correlated systems [27]. In our case, 
however, we recall that the shape of the PDF depends on the lag r and that for 
r 3> t s a Gaussian distribution is recovered. Several recent works have focussed on 
the variance of the fluctuations of the dynamics in glassy systems. In particular, 
Lacevic et al. [28] introduce in simulations of glass formers a time-dependent "order 
parameter" Q(t) that compares the system configuration at two times separated 
by t and then calculate the (normalized) variance of Q(t), Xii^)- These quantities 
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are analogous to our c/ and var(c/), respectively, since c/ measures the degree of 
overlap between sample configurations separated by a lag r. Indeed, X/t(t) is found 
to be peaked around the characteristic time of the final relaxation of Q{t), much 
as, in our experiments, var(c/) is peaked around t s . As a final remark, we point 
out that in most of the recent work on heterogeneous dynamics in glassy systems, a 
key role is played by spatial heterogeneity and spatial correlations between the re- 
arrangements. In our experiments, each CCD pixel is flooded by light issued by the 
whole scattering volume (see Fig. [Q. Accordingly any information on the spatial 
localization of the rearrangements is lost and only global temporal heterogeneity 
can be detected. We are currently working on a modified TRC scheme to overcome 
this limitation and achieve space- and time-resolved scattering measurements. 
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